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Validating Coherence Measurements Using 
Aligned and Unaligned Coherence Functions 


Jeffrey Hilton Miles 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


This paper describes a novel approach based on the use of coherence functions and 
statistical theory for sensor validation in a harsh environment. By the use of aligned and 
unaligned coherence functions and statistical theory one can test for sensor degradation, 
total sensor failure or changes in the signal. This advanced diagnostic approach and the 
novel data processing methodology discussed provides a single number that conveys this 
information. This number as calculated with standard statistical procedures for comparing 
the means of two distributions is compared with results obtained using Yuen’s robust statis- 
tical method to create confidence intervals. Examination of experimental data from Kulite 
pressure transducers mounted in a Pratt & Whitney PW4098 combustor using spectrum 
analysis methods on aligned and unaligned time histories has verified the effectiveness of 
the proposed method. All the procedures produce good results which demonstrates how 
robust the technique is. 


Nomenclature 


D delay time, sec 

d\ Yuen squared standard error of the aligned coherence values 

d 2 Yuen squared standard error of the unaligned coherence values 

f c upper frequency limit, f c = 1/2A t = r/2, Hz. (24000 Hz.) 

G xx {f),G xx (f) auto-spectra 

G xy (f) cross-spectrum 

j positive imaginary square root of — 

L y number of frequencies, / C /A/ = N/2 ( 2048 ) 
n number of coherence values used to compute A 
N1CORR corrected rotor speed, rpm 

rid number of disjoint (independent) segments used in spectra estimates of PW4098 test data,nd 
B e T to tai = 234 
random time series 
n s number of segments 

NP segment length, number of data points per segment (4096) 

r sample rate, samples/sec. (48000) 

s(t) time history 

t time, sec. 

Ttotai total record length, sec. (« 20 sec.) 

tyuen the 1 — a/2 quantile of Student’s t distribution with vyuen degrees of freedom 
Td(i) record length of segment i , sec. Td = = 4096/48000 = 0.08533 sec. 

X t \ sample trimmed mean of aligned coherence values 
X t 2 sample trimmed mean of unaligned coherence values 
X(f) Fourier transform of x(t) 

x{t) signal x time history 
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Y (/) Fourier transform of y(t) 

y(t) signal y time history 

Subscripts 

1.2 locations 1 and 2 

i running segment index 

Symbols 

a probability 

A / frquency step, l/T^Hz. (11.718) 

At sampling interval, 1/r (1/48000), sec. 

r] evidence dimension 

7 j? (/) estimated magnitude squared coherence (MSC) function 
Jxy(f) magnitude squared coherence function 
A confidence test parameter 

p, mean value 

i'Yuen Yuen number of degrees of freedom 

v degrees of freedom,^ = 2 (n — 1) 

a variance 

Superscripts 

* complex conjugate 


I. Introduction 

Sensor validation in the broadest sense is related to reliability. One must be able to determine that a sensor 
is or is not providing the correct signal. In addition, if the signal changes one would like a test that was able 
to distinguish between a sensor degradation from normal operating performance, total sensor failure, and 
changes in the signal due to changes in the signal source unrelated to the sensor. This document shows how 
conventional statistical methods 1-3 and a robust statistical method developed by Yuen 4-6 can be applied to 
identify Kulite sensor degradation or Kulite sensor failure, and separate these changes from changes in the 
signal . The question examined here is given two combustor Kulite pressure transducers how can the sensor 
signals be used to validate dynamic signal measurements. The sensor analysis scheme discussed herein was 
developed as part of a study of core noise from a Pratt & Whitney turbofan engine by Miles. 7 
The goals in analyzing the sensor data are: 

1. to minimize Type I errors, i.e. false alarms, and Type II errors, i.e. unable to detect an invalid signal 
when one occurs. 

2. to provide a rapid decision as to the existence of a faulty condition. 

3. to provide a test that does not depend on steady state (sample) statistics but rather statistics calculated 
from a time varying process. 

4. to provide a test that validates a pair of measurement devices. 

5. to provide a test that produces a single number that provides this information. 

Sensor validation plays an important role in automatic control and system monitoring within safety 
critical systems such as the aerospace and nuclear industry. 

A sensor validation scheme based on the auto-associative neural network is discussed by T.-H. Guo and J. 
Musgrave 8 for application to detect and recover from sensor faults in a simulation of the Space Shuttle Main 
Engine (SSME) control system. They reconstruct one or more lost sensor data using a redundant sensor set. 
Real time intelligent sensor validation using a Bayesian network for the detection of a fault in a set of sensors 
and a Bayesian network to isolate the faulty sensors discussed by Ibargitengoytia Ref. 9 
Another integrated sensor validation and fusion scheme is discussed by Wellington, Atkinson, and Sion Ref. 10 
One limitation of these methods is that validation of an individual measurement device is disregarded and 
system level validation is emphasized. 
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The following papers focus on individual measurement device validation through examination of patterns 
created by different sensor failures. 

Validation of thermocouple sensor signals based on their signal patterns has been discussed by Yung. 11 
Using the sensor signal from a dynamic process to validate a sensor has been discussed by Himmelblau and 
Bhalodia. 12 

A novel approach based upon using wavelet transforms to characterize the sensor signal in the time and 
frequency domain has been proposed by Ma, Zhang and Yan. 13 

A sensor fault detection and identification (FDI) scheme based on parameter estimation by means of an 
extended Kalman filter (EKF)is presented by Del Gobbo. 14 

Applying signal processing methods to signals from a single sensor for validation has been suggested by 
Amadi-Echendu. 15 

The paper is organized as follows. The theoretical background will be outlined next. Then the experi- 
mental results will be discussed. Finally, the conclusions are given. 


PW4098 Engine at C-ll Stand for EVNRC Phase 2 Tests 



Figure 1. Pratt &; Whitney test stand Cll, West Palm Beach Florida, with PW4098 engine and attached 
acoustic inflow control device also with and without aft acoustic barrier walls for EVNRC Phase 2 tests. 


II. Background 

The details of the aligned and unaligned coherence method are given by Miles. 7 Only an outline will be 
presented herein. The magnitude squared coherence (MSC) function, 7 xy (f) calculated using the complex 
cross-spectral density, G xy {f ) 1 and two auto-spectral densities at frequency f G xx {f) and G yy (f) is defined 
by: 
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(1) 


7 l y {f) = 


|G, y (/)| : 


G xx {f)Gyy{f ) 

The MSC function is discussed by Bendat and Piersol, 16-18 and Carter. 19 The MSC function identifies the 
degree of coupling between two time series. The normalization used insures the coherence has values between 
zero (the absence of correlation) and one (complete correlation). 

The MSC function estimation procedure used herein is based on signal processing algorithms developed 
by Stearns and David. 20 They refer to their spectral-estimation code as using a periodogram averaging 
method. Many hardware and software spectrum analysis systems use similar methods. The following set 
of equations will be used to discuss the calculation of the estimated MSC function by the periodogram 


averaging method. The estimated MSC function calculated using rid weighted Fourier transforms X{f) i and 
Y{f) i based on a segment record length Td is defined by 


w(t) 

x{t) 

y( T ) 

Mf) 

Yi(f ) 

7 UD 



w{t)x{t) 


w ( T )y( T ) 

rTd(i) 


/ x(t) exp (— j27r/r)c?r 

Jo 

/ y(r) exp (— j27r/r)dr 

Jo 

i £ r = i * i (/) fr (/) i 2 

Er= i E?=i ) 


(2) 


Here,* denotes the complex conjugate and the weighting function is w(r). 

As a calculated quantity the estimated MSC value is always nonzero. Consequently, nonzero values can 
appear even if the two series are independent. One needs a threshold value of MSC to decide if two series are 
related over some frequency range. Practical calculation of the estimated MSC show its statistical nature. 
Consequently, the statistical properties of the MSC have been studied theoretically by Carter 19 and using 
simulation method by Foster and Guinzy, 21 Benignus, 22 and Miles. 7 These studies indicate that if the true 
MSC is zero and the two time series, x n [t ) and y n (J) are independent then the probability, Pi, that the 
estimated MSC, |7x y | 2 , is in the interval, 0 < |7 xy | 2 < Ei is related to Ej by 


Ej = 1- (1 - P/) 1 ^"*- 1 ) (3) 

In addition, it can be shown that the expected value of the estimated MSC when the true MSC is zero, 
E^\lx.n.y n | 2 (/) | n s , |q| 2 = 0], is only a function of the number of independent blocks used to average the 
MSC. 


E[\% n yj 2 (f)\n a ,\ 1 \ 2 = 0\ = ±- (4) 

It'S 

The quantity |7x„ yn | 2 = Ei can serve as a threshold in deciding if two series are uncoupled. In addition, 
it is possible to calculate a confidence interval for |7x y | 2 (/) for any value of \'y X y\ 2 (f) using the conditional 
probability equations of Carter 19 and a method discussed by Wang et al. 23,24 

Unfortunately, the theoretical foundation for the statistical threshold value does not include the possibility 
that the samples used in estimating the MSC are overlapped. In addition, the use of different weighting 
functions, w(t), is common and might make the statistically derived threshold questionable. These type of 
questions lead to the use of simulation methods and surrogate methods to determine threshold values. 

The method of surrogate data analysis modifies a set of measurements to study its properties ( Faes, 25 
Thciler, 26 Schreiber 27 ). The method of surrogate data analysis has been used to determine a MSC threshold 
by Faes. 25 According to this method a set of surrogate series mimicking some properties of the original 
series but being otherwise uncoupled is generated. The threshold for zero coherence is then computed on the 
distribution of coherence estimates obtained from the surrogate series. For example, one of the procedures 
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used by Faes to generate a surrogate series was to randomly permute in temporal order the samples of the 
original series so that any temporal structure was destroyed in the surrogate generation. This procedure was 
devised for cardiovascular variability analysis where the coherence between the spontaneous fluctuations of 
heart period, systolic pressure, and respiration are used to quantify the amount of correlation between these 
oscillatory components and the reliability of transfer function estimates (Faes 25 ). 

This paper uses a novel measurement technique involving aligned and unaligned coherence functions and 
a novel data processing methodology to determine a threshold for zero coherence. First the measurement 
technique will be discussed. Then the data processing methodology will be outlined. 

A. Measurement scheme 

This aligned and unaligned coherence function measurement technique is described in great detail by Miles. 7 
The technique works because of time delay bias. Seybert and Hamilton 28 show that if one uses an average 
periodogram type calculation procedure to calculate MSC the time delay bias is 



where D is the delay time. 

This relationship has been used to explain time delay bias by Halvorsen and Bendat, 29 Seybert and Hamil- 
ton, 28 Bendat and Piersol 18 in chapter 11.2 (page 271) and Carter. 19,30 

The method of aligned and unaligned coherence is based on the following two characteristics of the 
average periodogram process when it is used to calculate the MSC : 

• When the time delay, D , exceeds the sample record length ,Td, the coherence of the random process, 
7xy(/)> i s the coherence of two independent random signals, 7 x nVn (f)- 

• If the signals contain tones using a time delay, D , greater than the sample record length ,T<j, does not 
change the coherence of the tones. 

B. Data processing scheme 

The unaligned coherence is identical to the coherence that would be measured if the signals from the two 
Kulites were independent. The unaligned coherence is used as a reference threshold for identifying the 
independence of two time histories. The coherence is studied without the frequency dimension using only 
the coherence magnitude. The magnitude of the coherence is a variable which has a distribution that 
expresses the difference between the aligned coherence and the unaligned reference coherence. In order to 
magnify the difference between the aligned and unaligned coherence a variable transformation is used to 
make the distributions more symmetric and normal. Various two parameter statistical hypothesis testing 
procedures using conventional statistical methods 1-3 and a robust statistical method developed by Yuen 4-6 
are applied to form a parameter that indicates sensor degradation, total sensor failure, or a change in the 
signal. The results are obtained since multivariate normal distributions play a pivotal role in the theory 
of multivariate statistical analysis which is used such diverse fields as control theory, stochastic systems, 
statistical signal analysis and multiple time series. 31-34 

Experimental data from Kulite pressure transducers mounted in a Pratt & Whitney PW4098 combustor 
is used to demonstrate the proposed method. This experiment is discussed first. Next the treatment of the 
Kulite signals is discussed. Then the data transformation and results of the statistical hypothesis testing 
procedures are discussed. 


III. Experiment 

To demonstrate the usefulness of the proposed aligned and unaligned coherence method, Kulite pressure 
measurements made in a Pratt & Whitney PW4098 combustor will be used. The measurements were made 
in a study of aircraft engine core noise conducted as part of the NASA Engine Validation of Noise Reduction 
Concepts (EVNRC) Program. One Kulite was at 127 degrees and the other was at 337 degrees. Kulite 
angles are measured clockwise from top dead center viewed from the rear or exhaust. The spectral estimate 
parameters are shown in table table 1 on the following page The signal processing algorithms used were 
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written in Fortran based on subprogram modules developed by Stearns and David 20 and were modified for 
this project so that the time series could be translated. In the calculations the segments were overlapped by 
50 percent. 

The test stand is shown in Fig. 1 
on page 3. The Kulites mounted in the 
combustor with water cooled jackets are 
shown in Fig. 2 on the following page. 

The three test cases shown herein are 
for test points with at 1622 rpm, 1999 
rpm, and 2600 rpm (N1 CORR). 

A. Combustor Kulites 

Selected plots of Kulite auto spectra and 
cross spectra are shown in Fig. 3 on page 
9. Corresponding aligned and unaligned 
coherence functions are shown in Fig. 4 
on page 11. 

Results of modal model analysis are 
discussed in Appendix B. Again, note 
that approaches to find a single number 
that describes the use of the aligned and 
unaligned coherence functions and sta- 
tistical theory to test for sensor degradation, total sensor failure or changes in the signal will be shown for 
a range of engine speeds. However, only data for three cases are presented 1622 rpm, 1999 rpm, and 2600 
rpm (N1 CORR). 

Given measurements from two sensors in a combustor to examine, we say a sensor is providing no 
information if the signals are uncorrelated with each other. In addition we suspect a sensor is providing an 
incorrect signal if the correlation changes abruptly in some irregular fashion. This paper uses the coherence 
function to determine if the signals from two sensors become uncorrelated or change in some strange fashion. 
The procedure used compares the normal measured coherence, called the aligned coherence with 

one calculated by time delaying one signal by an amount D greater than the sample interval Td . The value 
of D is chosen so that the two signals are not in the same sample interval, Td and appear to the processing 
procedure to be uncorrelated. In the paper a value of D approximately 1.5 times the sample interval is used. 
Thus to calculate the unaligned coherence we have 


Table 1. Spectral Estimate Parameters 


Parameter 

value 

Segment/block length i.e. Data points per segment, NP 

4096 

Sample rate, r, samples/second 

48,000 

Segment length, Td = NP/r, seconds 

0.08533 

Sampling interval, At = 1/r , seconds 

2.0833 x 10 -5 

Frequency step, A / = 1 /Td, Hz 

11.718 

Upper frequency limit, f c = l/2Af = r/2,Hz 

24000 

number of frequencies, Ly = / C /A/ = NP / 2 

2048 

Time delay, r = 6323/48000, seconds 

0.1317 

number of independent segments/blocks 

234 

overlap 

0.50 

Sample length, sec. 

20 


Tk(t) = Si(f)+7lfc(f) (6) 

rg(t — D) = S 2 {t — D) + m(t — D) 

The same signal processing algorithms are used for the unaligned time series. However, since the two series 
are now uncorrelated the ensemble average of the cross-spectrum is that of an uncorrelated process. 

The measured unaligned coherence, % r (/, D ) is not zero. For example, if tones are present in the time 
signal they will be in all ensembles and appear in the unaligned and aligned coherence. In addition to aiding 
the identification of persistent tones, the unaligned coherence provides a reference coherence. Small values 
of the aligned coherence greater than the unaligned reference coherence can be easily identified as significant 
since one knows what the uncorrelated coherence looks like. 

However, the major support it provides is in determining if a sensor is providing a useful signal. If the 
aligned and unaligned coherence are similar than the two sensors are uncorrelated and one sensor may have 
failed. For the experiment discussed herein, an examination of the two auto-spectra can be used to identify 
the sensor failure. 

The creation of the unaligned coherence has the same advantage as using a placebo in a medical study. 
In developing a drug or medical procedure due to factors such as the natural biological variability between 
individual patients and the placebo effect (i.e. a response due to therapy rather than the type of therapy) 
one cannot usually conclude that some therapy was beneficial on the basis of experience. Examples of 
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Figure 2. Kulite hardware mounted on PW4098 turbofan engine at Pratt & Whitney test stand Cll, West 
Palm Beach Florida for EVNRC Phase 2 tests. 


placebos are an injection of saline, a sugar pill, and surgically opening and closing without doing a surgical 
procedure. Creating the unaligned coherence from the same signals takes into account the variability due to 
the operating condition and provides a set of data one may compare with the aligned coherence to determine 
any difference. 

The data processing methodology will be discussed next and the statistical approaches used will be presented. 

IV. Transforming Data 

A number of schemes will be developed which produce a single number that provides knowledge about the 
degree of degradation from normal operating performance and the total failure of a pair of Kulite pressure 
transducers mounted in a combustor. This will show that having the unaligned and aligned coherence is a 
very robust procedure that produces information that can used by a variety of methods. 

The aligned coherence is a positive function due to broad band combustion noise, tones, and random 
noise. The unaligned coherence function is due to tones and random noise. One expects that at any 
particular frequency the magnitude of the unaligned coherence will be lower than the magnitude of the 
aligned coherence. 

When a coherence function is due to random noise, the coherence is low and the plot of coherence 
verses frequency consists of random points. However, the unaligned coherence can have a component due 
to tones in addition to one due to random noise. One result from the mathematical field of combinatorics 
is Ramsy’s Theorem which states that complete disorder is unlikely. No matter how jumbled and chaotic 
the arrangement the viewer creates a very highly organized and structured object within it. 35 Thus looking 
at a coherence plotted verses frequency due to random noise one sees patterns in the plot just as one can 
see shapes in the night sky created by randomly arranged stars. The presence of tones in the unaligned 
coherence make it even easier to see patterns when viewing the coherence as a function of frequency. 

In order to study the coherence without the frequency dimension and to include the whole range of 
coherence values measured in a frequency band just the magnitudes will be considered. The transformed 
magnitude of the coherence is posited as a dimension which has a distribution that expresses the difference 
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between the aligned and unaligned coherence function independent of frequency. In order to magnify the 
difference between the aligned and unaligned coherence a variable transformation will be used to make the 
distributions more symmetric and normal. The magnitude after transformation becomes an ’’evidence” di- 
mension. Larger values along this dimension above some threshold value provide evidence of correlation 
between two time series. Values below some threshold provide evidence that the two time series are inde- 
pendent Observations of the aligned and unaligned coherence have different distributions of the transformed 
function of the coherence because the ensemble average process shifts the aligned coherence along the evi- 
dence dimension to a higher level relative to the unaligned coherence. The evidence dimension used herein 
is 


Tj = 8 + 3.75 *log w (*il 2 ) (7) 

This scaling is selected to position the center of the unaligned coherence distribution near zero (77 = 0.0) since 
77 = 0 corresponds to 7 2 2 = 0.007356 = 10 -8 / 3 - 75 which is near the mean value of the unaligned coherence. 

The statistical treatment is based on the idea that the frequency label is more like an observation 
label. For each case the experimenter makes M observations i.e. the magnitude of the coherence and 
classifies them into a set of qualitative measurement classes. The results are independent of the order the 
observation are made in. Thus creating this statistical study the frequency associated with a coherence 
value is ignored. The study focuses on the statistical distribution of magnitudes of coherence values in a 
class. In order to get some idea of the frequency dependence this is done for five frequency intervals or 
classes: [0, 809], [0, 199], [211, 398], [410, 610], [621, 808] Hz. . The maximum number of points in the interval 
is limited since A / = 11.718Hz. This is a result of the use of a record length of 20 seconds and a selection 
of a reasonable number of ensemble averages. The number of points is reduced by neglecting points at the 
same frequency in the aligned and unaligned coherence which have levels which are closer than some set 
limit and might be due to tones. In the 800 Hz. interval there are typically sixty points and in the 200 Hz. 
interval there are typically sixteen points. The number of points seem sufficient to provide guidance as to 
sensor operation. 


V. Two Parameter Statistical Hypothesis Testing 

Many two parameter statistical hypothesis testing schemes are available . Selection of a scheme is based 
on knowledge about the statistical properties of the distributions such as if the population variance equals 
the sample variance. The following functions are calculated using standard statistical procedures based on 
the Student t test and Normal test 1-3 . The normal test is used if one may assume the distribution is normal 
and the population size is large. The Student t test is used if the population size is small and the probability 
curve is symmetric. 

They are used herein to compare the transformed aligned coherence magnitude and the transformed 
unaligned coherence magnitude distributions: 


A i ^aligned ^unalignedl 1 ’ 96 \/ cr aligned/ n + °unaligned/ n 

A 2 = I A 1 aligned “ ^unalignedl - tl ~ a / 2 ’ , '\/ (J aligned^ n + °unaligned/ n ® 

We call A a confidence test parameter. The Ai test is based on the normal distribution. The A 2 test is 
based on Student’s t distribution. The above formula applies if the number of samples n in the aligned and 
unaligned samples are equal. A more general t test is available if the two samples are not the same size . 3 
The quantity /i is the mean or expected value of the appropriate aligned or unaligned coherence functions 
over the desired frequency range. The quantity a 2 is the variance, which is the expected squared deviation 
from the mean of the appropriate aligned or unaligned coherence functions over the desired frequency range. 

Plots of these functions are shown in Fig. 5 on page 12. The curves were created using values of the 
parameters obtained by arithmetic averaging. 

When the function is positive the aligned coherence mean is not near the unaligned coherence and the 
sensor has not failed. However, in addition, the function provides information on signal changes and sensor 
degradation. Numerical values are given in table 1 and 2. 
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Figure 3. Auto-spectrum magnitude for Kulite (1) at 127 degrees and Kulite (2) at 337 degrees and cross- 
spectrum magnitude and phase between Kulite (1) and Kulite (2). 
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The third method uses a robust statistical method developed by Yuen 4-6 called the Yuen- Welch test 
This test examines the hypothesis that two independent groups have equal trimmed means and enables the 
computation of a confidence interval. This test is used if the population size is small and has outliers and 
the probability curve is nonnormal. 

The Yuen- Welch test states the 1 — a confidence interval for — Hn is 

A 3 = (Xtl(g) — Xt 2 (g)) ± ty uen \/ di + d ,2 (9) 

where X t i is the sample trimmed mean of aligned coherence values, X t 2 is the sample trimmed mean of 
unaligned coherence values, d\ is Yuen squared standard error of the aligned coherence values , ^2 is Yuen 
squared standard error of the unaligned coherence values, and ty uen is the 1 — a/2 quantile of Student’s t 
distribution with i>Y uen degrees of freedom. 

When A 3 is positive the hypothesis that the trimmed means of the aligned and unaligned coherence are 
equal is rejected. 

The Yuen procedure confidence interval is shown in figure 5 on page 12 and is marked as high (A]j~) and low 
(A^"). The Yuen procedure lower confidence interval value (A.j~) is in good agreement with the values of Ai 
and A 2 which are calculated using standard statistical procedures using the Student t test and Normal test. 

VI. Discussion 

Selected plots of auto spectra and cross spectra (figure 3 on the preceding page) and aligned and unaligned 
coherence (figure 4 on the next page) have been presented so that the sensor degradation, total sensor failure, 
and changes in the signal unrelated to the sensor operation as observed in the plots can be compared to the 
single parameter validity measurement discussed in the last section. These auto spectra and cross spectra 
plots when compared to the decision criteria proposed provide insight into the value of a particular approach. 
However, one might decide that using the approach presented in this section based on aligned and unaligned 
cross spectra provides sufficient information. Statistical theory adds no new information to the plots of the 
test cases. The statistical theory approach discussed does simplify and compress the results so that one 
figure can provide a test for sensor degradation, total sensor failure or changes in the signal. 

The results of Yuen’s confidence method and the standard statistical methods are shown in figure 5 on 
page 12. Results are shown for aligned and unaligned coherence function distributions of magnitude based 
on five frequency intervals : [0,809] Hz., [0,199] Hz., [211,398] Hz., [410,610] Hz., and [621,808] Hz. . 
All the procedures produce similar results. The fact that the robust statistical procedure of Yuen which 
compares trimmed means agrees with standard statistical procedures which compare mean values of aligned 
and unaligned coherence magnitudes shows that using the standard statistical procedures is feasible. 

The degradation of performance of the Kulite pair during the test is clearly shown in figure 5 on page 12 
The aligned and unaligned coherence are nearly the same at 2600 rpm (N1 CORR) where the confidence 
interval start to include zero and the values of Ai and A 2 go negative. Note that after failure the values 
go up. This is due to coherence between the signals from failed Kulite 2 and Kulite 1. A look at the Auto 
spectrum of Kulite 2 at 2600 rpm (N1 CORR) clearly show it has failed. 

The plane wave mode as identified as being in the 0 to 200 Hz frequency range Appendix B dominates 
at low engine speed (Appendix B). However, the most interesting point is that the plane wave combustion 
noise vanishes at higher rpm (N1 CORR) values while the noise in other modes remains. This is clearly 
shown by comparing the confidence interval plots in figure 5 on page 12. The curves in Fig. 5b covering the 
[0, 199] Hz. interval are below zero at 2000 rpm while the curves in Fig. 5c remain above zero until 2400 
rpm. 

The proposed tests go negative above 2500 rpm ( N1 CORR) as shown in Fig 5a when Kulite 2 fails as 
shown in figure 3 on the page before. Note however that the failure is not sudden but starts at 2250 rpm ( 
N1 CORR). 


VII. Conclusions 

A novel diagnostic procedure for sensor evaluation in a harsh environment based on aligned and unaligned 
coherence functions was compared with a method based on examination of auto-spectra. The new method 
is more effective in that it clearly indicates the signals are uncorrelated if a sensor fails. In addition, it does 
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Figure 4. Comparison of aligned and unaligned coherence between combustor Kulite (1) at 127 degrees and 
combustor Kulite (2) at 337 degrees. 
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this a priori .i.e without a prior model of the sensor. Furthermore, it enables one to have confidence that 
low values of coherence are significant. Examination of the coherence does not identify the failed unit. A 
failed unit can be identified using the Kulite auto-spectra. Examination of experimental data from Kulite 
pressure transducers mounted in a Pratt & Whitney PW4098 combustor by spectrum analysis of aligned 
and unaligned time histories has verified the effectiveness of the proposed method. 

Using a novel data processing methodology, statistical tests based on standard techniques and a robust 
statistical procedure by Yuen that uses trimmed means were examined to develop a single parameter measure 
of the validity of signals from the Kulite pair. All the procedures produce good results which demonstrates 
how robust the technique is. 

While the technique was developed to validate turbofan instrumentation, it might also be useful in 
medical fields where the coherence function is used to assess the significance of relationships of fluctuations 
in oscillatory components like heart period, systolic pressure, and respiration. 
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A. Auto Spectra and Coherence 


Figure 3 shows the magnitude of the auto-spectra and the cross-spectra and the wrapped and unwrapped 
phase angle of the cross-spectra for three engine speeds. 

At the 1622 rpm (N1 CORR)and 1999 rpm (N1 CORR) operating conditions the cross-spectra plots resemble 
those measured by Karchmer in a study of core noise from a YF102 turbofan engine Ref. 36,37 Karchmer 
showed that individual features of the combustion noise spectrum can be uniquely related to a specific 
propagation mode. It was found that the pressure field in the PW4098 engine combustor can be reconstructed 
in a similar manner that the pressure field in the YF102 combustor was reconstructed. 

The model used is based on the assumption that there are distinctive frequency bands over which the 
individual modes, including the plane wave can be treated as existing. Furthermore, for all modes but the 
plane wave mode the modes are assumed to overlap such that when mode m and mode m — 1 are present 
mode to — 2 is not present. 

Next we mention some of the results and illustrate them with a few examples. For the engine speed, 1622 
rpm (N1 CORR) case auto-spectra and cross-spectra are shown in Figure 3a and aligned and unaligned 
coherence functions are shown Figure 4a. Using the modal model results we see 

1. the plane wave mode A(0,0) below 187.5 Hz. 

2. the circumferential mode, A (1,0), between 187.5 Hz and 500 Hz. 

3. the circumferential mode, A (2,0), between 500 Hz and 800 Hz. 

4. the circumferential mode, A (3,0), above 800 Hz. 

At higher engine speeds similar model patterns appears. 

Examining the coherence plots in view of the results from the modal model suggest that in the 0 to 200 
Hz frequency band, the plane wave mode is very strong at 1622 rpm, 1703 rpm, 1750 rpm, and 1800 rpm. 
However, the plane wave mode starts to vanish at 1900 rpm and is gone at 1999 rpm. The presence of 
acoustic duct modes at higher frequencies as shown by the plots clearly indicates the sensors have not failed 
and that the noise source in the 0 to 200 Hz frequency range has changed. This is effect illustrated by Figure 
3b and Figure 4b where results for the 1999 rpm case are shown. 

The magnitude of the auto-spectra for Kulite (1) and (2) are similar until 2450 rpm. At that point Kulite 
(2) fails. At this engine speed the unaligned and aligned coherence become similar. In addition, the auto- 
spectrum G 22 does not resemble G\\. For the previous speed the two are similar. Above 2450 rpm the G 22 
spectrum is at least 20 dB higher than G\\. This behavior is illustrated in Figure 3c where auto-spectra and 
cross-spectra results for the 2600 rpm case are plotted. 

At each engine speed, the aligned coherence clearly shows the plane wave and circumferential propagation 
mode shapes. The unaligned coherence serves as a reference. It indicates, what the coherence would be if 
the Kulites were uncorrelated and the location of tones. The aligned and unaligned coherence are different 
until 2450 rpm. Figure 4 illustrates this behavior for the three engine speed cases discussed herein. A tone 
at 421.875 Hz is especially noticeable in many of the Figure 4 plots. At higher engine speeds above 2450 
rpm the aligned and unaligned coherence functions are similar indicating the two sensors are uncorrelated 
and one sensor may have failed. To identify the failed unit as Kulite (2), one must examine the auto-spectra. 
Fig 3c show a typical case. 
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Table 2. Normal distribution test based scheme. 

Ai = |/X S - Hn\ - Z\_ a !2\J &s/ n + a ll n 

95% confidence level a = 0.05 


rpm (N1 CORR) 

File 

[0,809] Hz. 

[0,199] Hz. 

[211,398] Hz. 

[410,610] Hz. 

[621,808] Hz. 

1622 

1756 

3.68 

2.79 

3.62 

2.70 

2.15 

1703 

1757 

4.02 

2.79 

3.70 

2.48 

3.12 

1750 

1758 

3.09 

1.25 

2.80 

1.78 

1.71 

1800 

1759 

3.01 

2.13 

2.34 

1.11 

.82 

1900 

1760 

3.57 

2.11 

3.79 

2.28 

2.49 

1999 

1761 

2.57 

-.66 

2.74 

1.81 

1.07 

2101 

1762 

2.90 

-1.19 

2.23 

-.55 

3.89 

2199 

1763 

2.88 

.15 

3.03 

-1.11 

3.14 

2304 

1764 

.37 

-2.05 

.38 

-2.04 

.07 

2400 

1765 

1.46 

-1.90 

2.51 

-2.19 

1.85 

2450 

1766 

1.90 

1.73 

.91 

-1.14 

.37 

2500 

1767 

.25 

-1.26 

-.92 

-1.92 

-.49 

2600 

1768 

-.28 

-.49 

-.96 

-1.82 

-1.05 

2700 

1769 

.99 

.07 

.53 

-.99 

-.64 

2743 

1770 

2.69 

3.70 

1.58 

-1.13 

.97 


Table 3. Student t distribution test based scheme. 

A 2 = | fls ~ Vn\ ~ tl-a/2^V a s/ n + a ?i/ n 

95% confidence level a = 0.05 


rpm(Nl CORR) 

File 

[0,809] Hz. 

[0,199] Hz. 

[211,398] Hz. 

[410,610] Hz. 

[621,808] Hz. 

1622 

1756 

3.67 

2.72 

3.56 

2.65 

2.01 

1703 

1757 

4.01 

2.72 

3.64 

2.40 

3.06 

1750 

1758 

3.08 

1.15 

2.70 

1.70 

1.63 

1800 

1759 

3.00 

2.06 

2.25 

1.02 

.64 

1900 

1760 

3.56 

2.04 

3.74 

2.09 

2.40 

1999 

1761 

2.56 

-.77 

2.67 

1.67 

.95 

2101 

1762 

2.88 

-1.27 

2.08 

-.67 

3.77 

2199 

1763 

2.87 

.03 

2.98 

-1.22 

3.07 

2304 

1764 

.35 

-2.26 

.33 

-2.15 

.01 

2400 

1765 

1.45 

-1.99 

2.43 

-2.30 

1.79 

2450 

1766 

1.89 

1.65 

.84 

-1.23 

.30 

2500 

1767 

.24 

-1.35 

-1.00 

-2.04 

-.56 

2600 

1768 

-.29 

-.55 

-1.05 

-1.91 

-1.13 

2700 

1769 

.98 

-.02 

.43 

-1.07 

-.73 

2743 

1770 

2.68 

3.64 

1.49 

-1.21 

.92 
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Table 4 Yuen procedure 
95% confidence level a = 0.05 


rpm(Nl CORR) 

File 

[0,80 

low 

9] Hz. 
high 

[0,19 

low 

9] Hz. 
high 

[211,31 

low 

38] Hz. 
high 

[410,61 

low 

L0] Hz. 
high 

[621, 8( 
low 

38] Hz. 
high 

1622 

1756 

3.81 

5.19 

2.77 

6.48 

3.12 

6.11 

2.52 

5.44 

3.35 

5.88 

1703 

1757 

4.22 

5.62 

2.81 

5.55 

3.42 

6.63 

2.96 

6.81 

2.88 

6.40 

1750 

1758 

3.46 

4.98 

2.05 

4.89 

3.72 

5.73 

1.96 

5.13 

1.09 

5.96 

1800 

1759 

3.67 

5.18 

1.92 

4.89 

2.34 

5.78 

2.25 

5.11 

1.83 

6.24 

1900 

1760 

3.50 

4.98 

2.32 

4.97 

3.16 

5.71 

1.64 

7.98 

1.79 

7.12 

1999 

1761 

2.70 

4.28 

-.30 

2.72 

2.51 

5.11 

2.78 

6.23 

1.56 

5.97 

2101 

1762 

2.93 

4.65 

-1.63 

2.08 

1.27 

7.58 

.46 

4.52 

3.18 

8.36 

2199 

1763 

3.03 

4.60 

-.14 

4.06 

3.01 

5.05 

-.40 

2.72 

2.93 

6.28 

2304 

1764 

.61 

2.07 

-1.64 

2.38 

.36 

2.41 

-1.06 

2.78 

.14 

2.68 

2400 

1765 

1.70 

3.53 

-2.06 

2.24 

2.58 

5.46 

-2.22 

2.14 

1.73 

4.23 

2450 

1766 

1.77 

3.52 

1.24 

4.62 

.36 

3.31 

-1.35 

2.14 

.10 

3.47 

2500 

1767 

.49 

2.07 

-1.53 

2.73 

-.62 

2.02 

-1.48 

2.66 

-.29 

2.59 

2600 

1768 

-.57 

1.51 

-1.15 

2.46 

-1.27 

3.29 

-1.50 

1.86 

-2.54 

1.25 

2700 

1769 

1.43 

2.86 

-.25 

4.33 

-.20 

3.74 

-.84 

2.26 

-.04 

2.70 

2743 

1770 

2.66 

4.14 

3.83 

6.97 

1.42 

4.67 

-1.23 

2.45 

.50 

3.29 
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